x1 = -1:0.01:1
y1 = x1.^2+1

x6 = -1:0.01:1
y6 = -x6.^2

x2 = -2:0.01:2
y2 = (x2+0.5).^2+1

x3 = -2:0.01:2
y3 = -(x3+0.5).^2

x4 = -2:0.01:2
y4 = (x4-0.5).^2+1

x5 = -2:0.01:2
y5 = -(x5-0.5).^2

hold on

%plot(x1,y1,'Linewidth',1,'Color',[0 0 0])
%plot(x6,y6,'Linewidth',1,'Color',[0 0 0])
plot(x2,y2,'Linewidth',1,'Color',[0.8 0 0])
plot(x3,y3,'Linewidth',1,'Color',[0 0 0.8])
plot(x4,y4,'Linewidth',1,'Color',[0.8 0 0])
plot(x5,y5,'Linewidth',1,'Color',[0 0 0.8])
xlim([-2, 2])
ylim([-1, 2])
xticks({})
yticks({})
xticklabels({})
yticklabels({})
%xlabel('k')
%ylabel('E')
%set(gca,'XMinorTick','on','YMinorTick','on')
%legend('paramagentic', 'up','down')
box off
hold off

set(gcf,'Position',[0, 0, 150, 150])